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Abstract 



■ A method of quantization of classical soliton cellular automata (QSCA) is put forward 
that provides a description of their time evolution operator by means of quantum circuits 

■ that involve quantum gates from which the associated Hamiltonian describing a quantum 
chain model is constructed. The intrinsic parallelism of QSCA, a phenomenon first known 

. from quantum computers, is also emphasized. 

O _ 

, Introduction. Soliton cellular automata (SCA), is a class of cellular automata|[l| operating 

on binary sequences with an updating rule function / for each cell, that depends on past and 
present time cells, the number of which determines the radius r of the automaton Contrary 
^ . to usual CA ( that evolve their cells by means of past time cells only ), SCA exhibit a variety 

I of evolution patterns |^, ^, ^ that is mainly known to characterize the temporal behavior of 

solutions of non linear PDE's||7|, namely: periodic evolution of particles (i.e localized groups 
of binary cells ), or solitonic type of scattering of digital particles, or even breathing modes 
of oscillations between particles. All these properties have motivated a number of suggestive 
^ ' applications for a new kind of computational architecture that will utilize these evolution pat- 

■ terns of SCA in order to provide a "gateless" implementation of logical operations!^. Towards 
a physical microscopic realization of these suggestions, envisaged in the context of the new 
paradigm of Quantum Computing^ it is plausible to formulate SCA in terms of Quantum Me- 
chanics and to investigate the possible quantum effects in their time evolution. To this end we 
put forward here the quantization of classical SCA, and point out their quantum parallelism. 

Quantization. Let the product Hilbert space n = ^i^jUx where Ui = span{\Q) 
Take n S N, and consider the subspace H C H where H = (i^''j^^_^?in^k, then define the vectors 
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and the dual vectors (a*| G i/, with orthogonormahty relation 



a* 1 6' 



n 



i=r \'^n-i\bn-i) ' (^nl^n) ' ^^=1 (o-n+j\^n+j) 



t+i ,t+i • ht ■ n 
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Recall that for {x, y) € the Kronecker delta function is defined as 5x^y = 1 ® x ® y, where 
e denotes XOR, i.e modulo-2 addition. Now take ''^^ = J.f^'^ \ 02'"+i, where O^r+i the 
2r + 1-fold null string, l]^^^"^ = 1.^^'^ \ {ao} where /({ao}) = O^^"*"^, the preimage of the 
null string, and define the one-to-one function[|5) / : Zg^*""^^ l}^^'^^ , as /(a*) = a*"^^, where 

where the updated bit takes the value a^"^^ = 1 ®\^^ a^'t^ ©j=o 

Introduce now the quantization of classical cellular automaton by means of the following 
quantization diagram. 



This diagram implies that: 
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The transition operator \J j implements in the space of qubits H the update rule / of the 
classical bits of the automaton, and acts trivially as the unit operator in the rest space. Since 
in order to establish the one-to-one property for / we have excluded the null string from 
its domain of values the associated operator f// € EndH, is a partial isometry in Hq, the 
orthogonal complement space of Hq = span{\0)}, assuming the decomposition H = Hq + Hq. 
To be more specific, let Uj = Z^^tg^*^"^^^ |/(^*) )(^*|' then its isometric property is due to the 
following relations: 

UfUf^ 
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a*eZr+\{a^} 



a*GZr+\{0} 
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where the preimage of the null string for e.g the simplest case of r = 2 is Og = 00100. 

Choosing a basis of vectors in H, e.g the computational basis B = {\xi) (8) • • • (X" |x2r+i) : 
G Zg''^^}, we can obtain a 2^*^+^ x 2^*^+^ matrix representation of Uf viz. 
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More explicitly if we partition the basis B into two orthogonal compliments corresponding to 
invariant and non invariant subspaces of C/j we will obtain for e.g the case of radius r = 2, the 
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transition matrix : 

Mu,) = ( ; ° ) . (6) 

where 1,0, stand for the square 16-dimensional unit and null matrices respectively, while o"^, is a 
N = 16-dimensional matrix with elements along the main antidiagonal {ai^^, a2^N-i, ■ ■ ■ , 07V,i) = 
(1, . . . , 1, 0), and all others been zero. 



Fig. 1. The factorization of Uj in terms of quantum gates for radius r = 2. 
Fig. 2. Quantum circuit implementation of FRT for three quantum BS's with r = 2. 

Hamiltonian model. In order to construct a Hamiltonian model associated with the QSCA 
that generates the total evolution of the automaton, and in view of the fact that the time step 
evolution operator is determined by unitary NOT viz. C/^ \i) = \1 (B i) and CONTROL — NOT(CN) 
viz. U^j^\ij) = \ii (B j) gate operators, cf. its explicit factorization in Fig. 1, we recall here 
first the Hermitian operators corresponding to those gates \ j) is abbreviated to The 

quantum negation of the i-th qubit is given by the matrix [/^ ~ ( 1 ^ ^ . ^'^^ trivially 

only in the i-th subspace, with associated Hermitian generator given by H^^ = + crj.), in 
terms of Pauli sigma matrices Q. For the conditional negation of CN gate with z-control and 
j-target qubits and i < j the unitary matrix C/^^ = |0 ) ( 0[ 1 + |1 ) ( 1[ (Ti, = exp{i-irH^(jj^) is 
generated by the Hamiltonian H^j^ = ^(1 — ct*)((T^ — 1). Correspondingly for i > j, C/^^y = 
IfSilO )(0|-|-(Ta;(8)|l )( 1| = exp(z7r//^^), and /f^^ = ^(1— (T^)((tJ, — 1). These gate operators can 
be utilized to factorize the U f that is used to update each qubit in a given string of qubits of a 
QSCA. To this end in Fig.l we construct the quantum circuit that provides such a factorization 
for the special case r = 2; the generalization to any r is straightforward. Specifically the input 
state vector in Fig. 1 reads |a*) = a*^_^2^ ('^n\^ 'Ai+i^ 'Ai+'i) ^ is tranformed by C// to 

I = rr |„A _ „t+l t+1 „t+l „t „t \ _ TjnTjn+2,njjn+l,njjn-l,njjn-2,n \ t\ ^>,p^p 

|o^+^) = 1 © © © © a^+i © 0^+2)- In Fig- 1 graphically the CN gate is indicated 
as a vertical line segment with a bullet (control qubit) and a crossed circle (target qubit) at 
its ends [10|. Then the total unitary evolution operator Uj = YVT ^/ = • • • ^i^i^o^ is the 
product of update operators for each qubit in the automaton starting with the first one Uq , 
from which all later ones are determined by shifting i.e ul_^_i = 1 © Uf . The total Hamiltonian 
Hf = lim(_+0 7"^ = J2i>o^i, is the sum of all Hamiltonians generationg each evolution 
operator {Uf , i > 0}. Since U/ = U\f 111=1 ^cn'^ Ilk=i ^cn'^' the corresponding Hamiltonian 
at each site i turns out to be Hf = H]^ + Y7k=i{-^c^N^ + ^c~n'^)- view of the generators of 
the quantum gates given above, the Hamiltonian model of the QSCA reads, 

i>0 k=l 

This can be interpreted as the Hamiltonian operator of a infinite quantum spin chain model 
with free boundaries and interactions ranged over r neighbors. Obvious modifications such that 
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QCSA with periodic boundary conditions, or higher dimensions are in the present formahsm 
appeahng and worth of further study. 

Quantum Fast Rule Theorem. We turn now to study QSCA in terms of the so cahed Fast 
Rule Theorem(FRT) that provides analytic tools for the classification and the study of types 
of dynamical evolution of classical SCA. In the limited space available here we shall confine 
ourselves to the special case of periodic particles and provide a quantum circuit of operators 
that governs the time development of periodic qubit-particles, i.e the quantum analogue of 
periodic particles of classical CA theory 

We recall that a particle is an integer multiple of r + 1 consecutive sites and it can be a 
collection of the so called basic strings (BS), which are r + 1 consecutive sites starting with a 
boxed site. If a particle consists of a single BS, we call it a simple particle. 
Then a basic theorem states the following^]: Let a sing le particle OA^A^ ■■■A^O with A^ 7^ 

0, A^ 7^ O, consisting of L basic strings A^,A^,---^ A^ (here O denotes the null BS consisting 
of r + 1 zeros), and let ® denotes sitewise XOR among bits of the BS. Suppose that the BS's 
A\ © A^, A^ ®A^,--- , A^~^ © A^, A^ , contain IqJi, - ■ ■ ,Il I's respectively. Then if such a 
particle does not split or loose any BS's from its right end at all times during its evolution, then 
at times t = I0 + I1 + . . .Im, m<k<L, it becomes A"^+'^ © (yl'^+2^m+3 . . . ^i^o^i . . . j^my 
Especially for k = L, i.e at time t = Iq + li + . . .Il it returns to its initial state A^ A? ■ ■ ■ A^ . 

To illustrate the workings of the quantum FRT i.e the quantum adaptation of the preceding 
theorem, let us consider an initial particle of L = 3 classical BS's A^ A^' . Such a classical 
particle by means of the quantization map p becomes the qubit word \ A^A^A^OOOO). Let us 
assume that the classical word labelling this state vector evolves under the previously stated 
conditions so that the classical FRT applies and assures that it has a periodic behaviour. 
Then the quantization map induces into the quantum state vector the periodic evolution of 
its classical label word. Namely the quantum state vector evolves periodically, and we seek 
to determine the evolution operator that implemends this periodic motion. Assuming we 
have a sufficient number of null quantum BS's from each side of the initial particle state 

1. e \A^A^A^OOOO), let us consider the following two operators. First, the operator Vq = 
^xe'Z.*^^^ jO'""'"^ ){x\, which projects each r + 1-fold tensor product BS state vector onto the 
r + 1-fold null state vector. With the upper index / denotes the fact that the (r + l)-fold 
input vector is placed in the I-th position of a given chain of tensor product of state vector, 
and that Vq acts non trivial only in that I-th subspace. Second, let us define the collective CN 
gate operator U^fj^ = Uq^^ '■+1 . . . Ujj^'^Ujj^^ ^ with control BS placed in the I-th position and 
target BS placed in the J-th position. This collective CN gate is actually defined as a succesion 
of (r -|- 1) CN gates for the corresponding qubits of the (r -|- l)-fold tensor product control and 
the target BS states vectors. 

In Fig. 2 we indicate graphically the sequence of operations that tranform the initial 
state vector \A^A^A^OOOO) into the final state \OOOOA^A^A^), indicating in this way its 
periodic propagation. Each parallel wire in the figure represents an (r + l)-fold tensor product 
of qubits and the operators Vq are indicated by boxed zeros, while operators U^"^^ are indicated 
a double circled CN gate symbol. Explicitly the sequence of operations of the quantum circuit 
in Fig. 2 reads: 

A^A^A^OOOO) 
OA^(S)A^A^®A^A^OOO) = Vlu}.%U^\U^ 



A^A^A^OOOO^ 
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ooA^ e A^A^A^ e A^ooj = vlul\ul%ul\ OA^ e A^A^ e A^A^OOO^^ 
oooA^A^ e A^A^ e A^o) = v^u'^%u^%u^\ looA^ e a^A^a^ e .4200 J 
oooo^^^2^^) = r^u^^N^^%u^% oooA^A^ e a^a^ ^ ^s^v ^ 



Quantum parallelism. The quantum description of the cehs of a C A in terms of qubits allows 
for having a superposition a |0) + 6 |1) at each cell. According to the standard interpretation of 
Quantum Mechanics this combination means that we have the state |0)(|1)) with probability 
|ap([6p). This is the inherent probabilistic character of QCA that entails two important 
advantages of QCA over classical probabilistic (noisy) CA (see e.g [||]). First, there is an 
exponential overhead in storing cell values in QCA over classical CA: an r-radious probabilistic 
SCA needs at each discrete time step to store = 2^'^"'"^ input words in order to process them 
later on and to update the current cell value. On the contrary a QSCA may form and admit 
as input a superposition state vector lijj) = '^12^^-^*'^^+'^ \^) = '^n^xi^'Z.'^''^^ ' ' '^Sr+i)) 
which with resources linear in r i.e 2r + 1 qubits only can constructs an exponential in r 
register of 2^''"'"^ equiprobable input states. Second, there is an exponential overhead in the 
updating time of the cells. Once a superposition of all input states has been prepared as 
the single state vector 1-0), then by acting the linear evolution operator Uf only once on it 
i.e Uf\ij)) = E^jg^'^^+i ^/ 1^) = ^xeZ*^"^^^ l'^^^-^^' ^® ^^'^ update all the state 
vectors simultaneously. This same exponential acceleration in storing space and processing 
time of quantum information has been known in quantum computation, where it is utilized 
in a number of tasks such as fast quantum algorithms, information encoding, cryptography 
etc§. 

Harnessing the space-time resources made available by the quantum parallelism of QCSA 
in various applications is an open challenge. The research program initiated here will require 
further studies of QSCA both as a new computational machine and as an implementable 



physical system. Some of these questions however will be addressed elsewhere|ll]. 

Discussion. We close by discussing some prospects of our work. Towards a realization of 
QSCA as a physical process that could be implemented experimentally and so would potentially 
provide a novel quantum computational machine, we may construct an optical circuit made 
of elementary passive optical elements i.e light beam splitters and phase shifters that are 
assembled so that they implement the operator C/j of a QSCA. For that purpose we may utilize 



the theorem of embedding the SU{2) Lie group into SU{N), which recently|12] has been used 
in the form of expressing every unitary matrix by a sequence of embedded 2x2 elementary 
unitary matrices of only two types : one matrix that realizes a beam splitter with transmission 
and reflection coefficients determined by the matrix elements and one that similarly realizes 
a phase shifter of a classical propagating light wave. By virtue of this embedding we may 
factorize the Uf matrix of a r = 2, say QSCA into elementary matrices of beam splitters and 
phase shifters. Assuming a two-state encoding of the qubits of the QSCA cells (e.g taking the 
vertical and horizontal polarization states of a laser beam as |0 > and |1 >), we may construct 
the quantum optical analog of the Uf of the QSCA. 

Finally, we should mention a number of ramifications of our QSCA formalism that are 
currently also under investigation, namely: QSCA beyond the binary case ( this would involve 
higher dimensional representations of the discrete Heisenberg group ) ; differential versus matrix 
formulation of QSCA (i.e partial differential operators acting on a space of multivariable 
complex polynomials realizing the respective action of the Uf matrix ); and also the case of 
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noisy QSCA (i.e QSCA evolution that occurs in the presence of of quantum mechanical noise, 
exemplified by errors due to random bit-flipping and phase shifting in the qubits of the QSCA 
string ); this would require the description of QSCA in terms of trace-preserving operators 
instead of the unitary operators as in the noiseless case, and the respective quantization of 
the classical SCA to be carried out not by Hilbert-space state vectors but by employing the 
/9-density operator formalism of Quantum Mechanics. 
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